Fast quantum noise in Landau-Zener transition 
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We show by direct calculation starting from a microscopic model that the two-state system with 
time- dependent energy levels in the presence of fast quantum noise obeys the master equation. The 
solution of master equation is found analytically and analyzed in a broad range of parameters. The 
fast transverse noise affects the transition probability during much longer time (the accumulation 
time) than the longitudinal one. The action of the fast longitudinal noise is restricted by the shorter 
Landau-Zener time, the same as in the regular Landau-Zener process. The large ratio of time scales 
allows solving the Landau-Zener problem with longitudinal noise only, then solving the same problem 
with the transverse noise only and matching the two solutions. The correlation of the longitudinal 
and transverse noise renormalizes the Landau-Zener transition matrix element and can strongly 
enhance the survival probability, whereas the transverse noise always reduces it. Both longitudinal 
and transverse noise reduce the coherence. The decoherence time is inverse proportional to the noise 
intensity. If the noise is fast, its intensity at which the multi-quantum processes become essential 
corresponds to a deeply adiabatic regime. We briefly discuss possible applications of the general 
theory to the spin relaxation of molecular magnets. 

PACS numbers: 03.65.-w;33.80.Be;32.80.-t; 
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I. INTRODUCTION 

Landau-Zener (LZ) theory^,— ^ deals with a two-state 
system, whose levels vary in time by sweeping of some 
physical parameter and would cross at some moment. 
Transitions between the two states are induced by the 
transition matrix element A, which in static or adiabatic 
case causes the level repulsion (avoided crossing accord- 
ing to the Wigner-Neumann theorem). The LZ theory 
is one of basic dynamic problems of quantum mechanics. 
It has a broad range of applications from the collision 
theory and quantum chemistry to disordered solids and 
molecular magnets. 

The LZ transitions in noisy environment recently at- 
tracted new attention in connection with the problem 
of the qubit decoherence^)^. New experimental realiza- 
tion of qubits- rose again the question to what extent it 
is possible to minimize the decoherence simultaneously 
maintaining sufficiently strong coupling to the external 
signal. 

The problem of the LZ transition in noisy environ- 
ment has relatively long history. One of the first con- 
siderations based on ideas of stochastic trajectories be- 
longs to Kasunokii. In the pioneering work^ Kayanuma 
calculated the transition amplitude in the presence of a 
fast transverse Gaussian classical noise with a specific 
(exponential) two-time correlation functional. This so- 
lution was simplified and extended to general shape of 
correlation function by one of the authors (VP) and N. 
Sinitsyn-il. In the same work a situation in which the 
transitions are produced by noise as well as by regular 
Hamiltonian was considered. Pokrovsky and ScheidU^ 
calculated the two-time correlation function of the tran- 
sition probabilities for the LZ system subject to a fast 
classical transverse noise. Longitudinal noise was consid- 
ered by Kayanumai^i, who proved that strong fast longi- 



tudinal noise does not change the LZ transition proba- 
bility. Gefen et al}- and Ao and Rammer^^ considered 
more wide range of parameters and found the situations 
in which the noise changes the transition probability. In 
the workii a rather detailed analysis of different limiting 
cases of temperature, coupling to the phonon bath, its 
spectral width and sweeping rate was presented. There 
occurred a controversy between the worksi^ ancU^. Gen- 
erally, there is no complete agreement between different 
authors on what happens in the adiabatic regime (very 
slow sweeping) in the presence of the longitudinal noise. 
Motivated by this disagreement Kayanuma and Naka- 
mura performed a comprehensive analytical and numeri- 
cal study of the LZ transition in the presence of longitudi- 
nal nois&i^. In particular they obtained a formula for the 
case of strong decoherence which is valid in both low- 
temperature and high-temperature limits. In all these 
works the quantum nature of the longitudinal noise was 
taken in account. 

Despite of significant progress a complete theory of 
the LZ transition in noisy environment still does not ex- 
ist. Theoretical works considered either quantum longi- 
tudinal noise with transitions originated from the reg- 
ular transition matrix element or the classical trans- 
verse noise. Quite recently Wubs et alJ^ have found 
a beautiful exact formula for the transition probability 
of the 2-state system interacting with the phonon bath 
at zero temperatureii. The noise had both longitudinal 
and transverse components. Their correlation and quan- 
tum nature were substantial. No limitations to the noise 
strength and spectral width were assumed. However, the 
limitation of zero temperature (phonon bath is in the 
ground state) does not allow to extend these results to 
more realistic situations. 

The purpose of this article is to present a theoretical 
description of the LZ system subject to a fast quantum 
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noise which has both transverse and longitudinal com- 
ponents. It is not yet complete theory, since it does not 
cover slow and intermediate noise, but in its range of 
applicability it allows to understand clearly all relevant 
physical regimes and phenomena. We will show that, due 
to the fastness of the noise, the LZ transition in the pres- 
ence of the longitudinal noise and the transitions due to 
the transverse noise are separated in time, whereas the 
correlation between the transverse and longitudinal noise 
leads to a renormalization of the regular transition ma- 
trix element in the LZ Hamiltonian. For a moderately 
strong transverse noise we derive master equations gov- 
erning the population of the two states and study their 
solution. If the transverse noise is strong and also fast, 
the 2-state system falls into adiabatic regime. The pop- 
ulation of levels comes to the equilibrium with the spin 
bath if the bath is in the state of thermal equilibrium. We 
argue that a very strong noise is classical and adiabatic. 
In this situation, as it was shown in^^, the populations 
of the two levels are equal. 

The plan of the article is as follows. In section 2 we 
introduce the Hamiltonian and characterize the noise. In 
section 3 we present simple heuristic arguments resulting 
in master equations. In section 4 we derive the master 
equations starting from microscopic Hamiltonian for the 
case of the transverse noise only and zero LZ transition 
matrix element. In section 5 we derive the renormaliza- 
tion of the regular transition matrix element due to cor- 
relation of longitudinal and transverse noise. In section 
6 we analyze the influence of the longitudinal noise. In 
section 7 find the solution of master equation and study 
it. In section 8 we match the solution of the master 
equation with the solution of the LZ problem without 
transverse noise. Section 9 contains discussion and con- 
clusions. Here we compare our theory with that by Wubs 
et ali^. We briefly analyze possible applications of our 
theory to molecular magnets. 

II. STATEMENT OF THE PROBLEM 

We consider a 2-state system interacting with a noisy 
environment. The latter is a large system (bath) with 
a stationary density matrix. We neglect the influence 
of the LZ transitions onto the state of the bath. For a 
definiteness we will speak about the phonon bath, though 
it can include other Boson excitations like spin waves, 
excitons, photons. Then the total Hamiltonian of the 
system can be represented as follows: 

H = H2 + Hb + H,„t (1) 

The term H2 in equation ([1]) represents the two-state 
system: 

H2 = -^a^ + Aa,, (2) 

where and ct^ are Pauli matrices and Q (t) is the 
time-dependent frequency or the energy difference be- 



tween the so-called diabatic levels turning into zero at 
t = 0; A is the regular transition matrix element. If t 
would be not time but some parameter of the Hamilto- 
nian, then non-zero A provides repulsion of the adiabatic 
levels (the Wigner-Neumann theorem on avoided levels 
crossing). For brevity we will call further the regular 
transition matrix element A the LZ gap. Usually the lin- 
ear approximation for the frequency 17 (t) = fit proposed 
by Landau and Zener is acceptable, but sometimes it is 
necessary to go beyond this approximation. Namely, in 
real experiment the sweeping of f2 (t) stops at some fi- 
nite value, which can be not large in the frequency scale 
of the problem. Therefore we will keep notation (t) 
throughout this article. The term Hi, in equation ([1]) is 
the phonon bath Hamiltonian: 

i7b = ^C.q6^6q, (3) 

q 

where 6q and 6^ are the operators of the phonon annihi- 
lation and creation; ujq are the phonon frequencies and q 
are their momenta. The interaction Hamiltonian reads: 

Hint = u\\(T^ + uj_a^ (4) 

The Hermitian operators U|| and u± responsible for the 
longitudinal and transverse noise depend linearly on the 
phonon operators. Each of them is a sum of two Hermi- 
tian conjugated operators containing either phonon an- 
nihilation or creation operators only: 

Ua ^Va + vi; Va = -j^^ 9a (q) bq, ||,-L (5) 

q 

where (q) are complex coupling amplitudes and V is 
the volume of the system supporting phonons. Quan- 
tum character of the noise manifests itself in non- 
commutativity of operators rj^ and 77 J,. The problem 
consists in calculation of transition and surviving prob- 
abilities for the 2-state system at a fixed noise density 
matrix. In the absence of the noise, the transitions am- 
plitudes constitute the LZ transition matrix belonging to 
the SU(2) group and depending on the dimensionless LZ 
parameter ^ A^/fJ: 

nz - ( ) (6) 

V2^exp(-^+zf) 

VrEiT i-nLz) 

If 72^2 is small, the system with the probability close 
to 1 remains in initial diabatic state; if 7^^ is large the 
system with probability close to 1 proceeds along the 
adiabatic state, i.e. changes the initial diabatic state to 
the alternative one. A characteristic time necessary for 

the LZ transition is t^z = max ( A/f2, {}~^^^ ] . 
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The influence of the noise onto the 2-state system is 
completely described by the noise correlation functions: 

(t) vl it')) ^^T.9a (q) 9*0 (q) ("q + 1) ^-''(''-O 



(7) 



(vi (i) Vp it')) ^lpY.9a (q) 9*p (q) n^e^^^i'-'') (8) 



Here riq = (&]j6q) are the average phonon occupation 
numbers and (...) means averaging over the phonon bath 
ensemble, not necessarily in thermal equilibrium. The 
Fourier-components of the correlation functions (|7l8p 
read: 



\ 2,71 

^/s) ^irY^da (q) g; (q) ("q + l)S{uj- wq) (9) 



V 



(viVp)^ ^y^9a (q) 9p (q) ^^q'^ + ^q) (10) 

q 

Note that one of the two correlators contains only posi- 
tive, whereas the second one contains only negative fre- 
quencies. If the noise is in equilibrium at temperature 
T, the Fourier-transforms of correlation functions obey a 
simple relation (uj > 0): 



(11) 



Let us denote LUg the range of frequencies in which 
the coupling coefficients ga (q) do not vanish. If the oc- 
cupation numbers Hq are of the same order of magni- 
tude for all states within this region of frequencies, then 
ujg determines the spectral width of the noise Alo. In 
some cases, for example at low temperature T <ti Wg, 
there appears a second, smaller scale of frequency (T). 
The noise correlation time is r„ = 1/Auj. It is differ- 
ent for the two correlation functions (rj^ (t) rj^^ (t')^ and 

{vti (*) Vf} {t')) at low temperature and it is equal to Ljg^ 
for both at high temperature. By definition the noise is 

fast if T„ < min (Cl'^^'^ , A'^^ or Aw > max (n'^/'^.A^. 

Besides its spectral characteristics the noise is character- 
ized by its strength. The most natural measure of the 
noise strength is the average square of its amplitude: 



7 J2 (q) 



(2n, 



(12) 



the perturbation theory is generally invalid at accumu- 
lation time scale determined in the next section, we will 
show that it works during sufficiently small intervals of 
time still longer than r„. Most of our results relate to 
the moderately strong noise. The noise is called strong if 
i^a (t)) ^ '''n^- During another time scale, which wc call 
the relaxation time t^q — {(u\ (t)) r„) ^ (seei^, where it 
is called phase relaxation time)the probability of finding 
the system in state 1 or 2 changes significantly. When the 
noise is strong, the relaxation time becomes less than the 
noise correlation time t„. It is convenient to introduce 
dimensionless coupling function Acq — \ga (q)| / (w^a'^), 
where a is the lattice constant. In condensed matter 
systems the values Aaq never become large. Migdal— 
argued that large coupling constants would lead to the 
lattice instability and reconstruction. Though his argu- 
ments related directly to the electron-phonon coupling, 
his idea is very general. If Aaq are not large, a large 
value of (m^ [t)) can be reached only if phonon oc- 
cupation numbers riq become large. In the equilibrium 
case it means that the temperature must be large. Very 
strong noise is classical irrespectively of its specific den- 
sity matrix. An analogue of the LZ parameter for the 
noise reads": 7„„ = (u^ {t)) /VL. If it is small, the noise 
brings only a small perturbation to the LZ picture; if 
it is large, the occupation numbers of the 2-state system 
follow adiabatically the instantaneous value of frequency. 
We will return to this point later. 



III. HEURISTIC APPROACH 

We start with an auxiliary problem specifying A = 
and r/|| = 0, so that all transitions are only due to the 

transverse noise. Since r„ <C the instantaneous 

frequency of the 2-state system does not change during 
correlation time of the noise and it can be considered as 
a constant. Thus, it is possible to calculate the instan- 
taneous rate of transition probability using the standard 
quantum mechanical technique for transitions between 
stationary energy levels. 



FIG. 1: Feynman graph for a 3-phonon process. Thin solid 
lines correspond to the state 1; thick solid lines correspond to 
the state 2; dashed blue lines correspond to phonons. 



The noise is weak if (u^ (i)) ^ f2. Weak noise can be 
accounted as a small perturbation to the LZ result. We 
call the noise moderately strong if it obeys the inequal- 
ity: (u^ (<)) <C T~^. Though for moderately strong noise 



For moderately strong noise the transition probability 
can be calculated in the first Born approximation. In- 
deed, according to the Fermi golden rule the rate of tran- 
sition probability per unit time before the level crossing 
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reads: 



Pi^2 (t) = 27r<^?7^r/^^ 



(13) 



In the framework of the considered model the next after 
single-phonon is three-phonon transition shown in Fig. 
1. Its contribution to the transition probabihty reads: 



pfi2 (t) = 27r 



0J1UJ2 



X herm. conj. 



(14) 



The 3-phonon contribution can be neglected if the noise is 
moderately strong. In the same approximation it is pos- 
sible to neglect the correction to the transition frequency 
due to the interaction with the phonon bath. Thus, the 
occupation numbers of the diabatic states Ni^2 at nega- 
tive time obey following master equation: 



-n(t)l \ /|o(t)| 

(15) 

which must be complemented by conservation law A^i -|- 
N2 — 1- For positive time equation ([T3]) must be modified 
as follows: 



iVi = 27r - 



\n(t)\ 



|-o(t)| 



(16) 

The noise produces transitions as long as its spectral 
width exceeds the instantaneous frequency |ri (t) |. The 
accumulation time estimated from this requirement is 

Tacc = (^^tn^ ■ Sincc the noise is assumed to be fast 
the accumulation time Tacc is much longer than the noise 
correlation time r„. In real experiment the sweeping of 
frequency may stop or saturate before the accumulation 
time is reached. The master equations enable one calcu- 
lating the occupation numbers at any time rather than 
asymptotically at t — s- cx). The accumulation time is 
also much longer than the LZ time t^z- Therefore, it 
is possible to neglect the action of the noise during the 
LZ time interval {—tlz,tlz) and neglect the LZ gap A 
beyond this time interval. It means that the action of 
the fast transverse noise and of the LZ gap are separated 
in time as it was earlier shown for classical nois&ii. The 
solution of the LZ problem without transverse noise and 
the noise transition problem with zero LZ gap A should 
be matched at some intermediate time. It will be done 
in Section Eini 

The action of the fast longitudinal noise is very dif- 
ferent from that of the transverse one. The longitudi- 
nal noise does not produce transitions in the absence of 
the LZ gap. Therefore, its action is effectively limited 
to the LZ time interval. In the next section we demon- 
strate that the fast longitudinal noise must be sufficiently 
strong to produce a substantial change in the LZ tran- 
sition probability. Namely it must satisfy an inequality 



u|y > Q/ (Ar„) '3> fJ. An analogous criterion for the 

transverse noise is much more liberal: (u^) > fl. For a 
comprehensive analysis of the longitudinal noise action 
we refer the reader to the cited articlesi^r^ji^. Beyond 
the LZ time interval the classical longitudinal noise mod- 
ulates the transverse noise by a factor exp (^—i J^^ U||C?t^ . 

Correlation functions (^V^l (t) {t')'^ must be substi- 
tuted by (i) 77 (f) exp l^—i Jj* U|| c?T^ y Neglecting 

the correlation between longitudinal and transverse noise 
and employing the Gaussian statistics of the noise, one 
can express the latter correlator as follows: 



V± it)v{ (Ocxp (^-ij^^ n^r^) = {v± it)v{ it')) 
dtidt2 (ii)uj| {t2)) 



X exp 



(17) 



The transverse noise correlator decays rapidly when 
the modulus of difference \t — t'\ exceeds r„. Therefore, 
the value in the exponent in the right-hand side of equa- 



Tt < 1. This esti- 



tion (fT7|) can be estimated as 

mates shows that the longitudinal noise can be neglected 
beyond the LZ time interval. In the next section we con- 
sider this question in more details. 



IV. DERIVATION OF MASTER EQUATIONS. 



to 



to 



FIG. 2: An example of a term in the perturbation theory. 
Points correspond to vertexes Vi{tj). 

Our goal is to find the dependence of the occupation 
numbers Na (a = 1,2) on time. The same problem 
can be formulated as calculation of the average value 
of projectors Pa = \a) {a\ = ^ {1 ± az)^r^- We con- 
sider the case A = 0, W|| = 0. The calculation will be 
performed in the interaction representation with the di- 
agonal time-dependent unperturbed Hamiltonian Hq = 
-^a, + Hi,^ (|1) (1| - |2) (2|) + Hh and the in- 

teraction Hamiltonian V = u^a^ = uj_ (|1) (2| + |2) (1|). 
Being transformed to the interaction representation, the 
interaction Hamiltonian depends on time as follows: 

Vi it) = UA_ (t) (|1) (2| e-'^4 f^(^)'*^ + |2) (1| "W-*^) , 

(18) 
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where 

(t) = e'^o'-^-'^'^u^ (to) e-'-f^'(*-*«) (19) 

The time-dependent occupation numbers can be ex- 
pressed in terms of the evolution operator [// (t, to) in 
the interaction representation: 

TV, (t) = Tr {t, to) PcUi {t, to)] (20) 



\ * \ — " < — 1^— 

I 

^— 



FIG. 4: Elementary graphs, a) b) without phonon line cross- 
ing; c) d) with phonon line crossing. 



where Pq is the initial density matrix which is the direct 
product of two independent density matrices Pq = P2Pby 
where the first factor is the density matrix of the 2-state 
system and the second once is the same for the bath. 
For calculation or these averages we employ a simplified 
version of the Keldysh-Schwinger technique^,— used 
already for a similar purpose iii^^^. Each of the 
two evolution operators is presented as a series of time 
ordered integrals. A general term of such an expansion 
contains a product of two multiple time integrals. With 
each time variable tk & vertex Vi (tk) is associated. The 
product of vertices is ordered chronologically in Ui {t, to) 
and antichronologically in Uf^ (t^to). All operators 
of Vi (tk) belonging to Uf^ {t,to) are located on the 
left ("later") than corresponding operators belonging 
to Ui(t,to). A particular contribution is graphically 
depicted in Fig. 2 before the averaging over phonon bath 
is performed. It consists of two lines both starting at 
^0 and ending at t. The upper line symbolizes Ui (t,to) 
and the lower one symbolizes U]~^ {t,to)- Vertices on 
these lines correspond to the operators Vj (tj,). Each 
vertex contains one phonon operator u{tk) and changes 
the state of the 2-state system. The presence of the 
projection operator Pa in equation (|20p implies that the 
state closest to the final time t on both lines must be \a). 
Applying the Wick's rule for phonons, one should form 
all possible pairing of (phonon) noise lines^^i. However, 
in the case of fast moderately strong noise the main 
contribution to the occupation numbers comes from the 
graphs without crossing or overlapping of the phonon 
lines. An example of a graph without overlapping is 
shown in Fig. 3. In comparison with elementary graphs 
without phonon line overlapping or crossing (Fig. 4a, b) 
the contributions of the graphs containing overlapping 
or crossing (Fig. 4c, d,) have additional small factors 
of the order of (uj_)t^ and can be neglected if the 
transverse noise is moderately strong. Indeed the time 
interval between the ends of each phonon lines is about 

«> I > ■ 

\ ! 
\ I 
\ I 
\ I 

^ — ( — 1^ — h ^ — t — — 

FIG. 3: A typical graph without phonon line crossings domi- 
nantly contributing to the survival and transition probability. 



To derive the master equation we consider a set of 
graphs differing the occupation numbers Na {t + At) 
from those for Na (t) with At satisfying the following 

strong inequahty t„ <C At ^ ((^1) '''") ■ First, there 
are graphs with phonon lines connecting the interval 
{t,t + At) with the interval {to,t). Their contribution 
can be neglected since it is relatively proportional to a 
small ratio r„/At. The contribution of k non-overlapping 
or intersecting noise lines inside the interval (t, t + At) is 

proportional to r^At)*^ <C 1. Therefore, the dom- 

inant contribution to the set comes from graphs contain- 
ing exactly zero or one line inside the interval (t, t + At). 
Note that it is a coarse-grain description: the master 
equation is invalid at time scale t„ and shorter. The 
graphical equation connecting Na {t + At) with Na (t) is 
shown in Fig. 5. 




FIG. 5: Graphic equation connecting A'^^ {t + At) and Na (t) 

To find their analytical expression, we calculate the con- 
tribution of the 3 elementary subgraphs for a = 1 shown 
in Fig. 6. 



^ B t2 tl 

FIG. 6: Three elementary graphs with one phonon line in the 
interval t,t + At. 

They read: 

t+At t+At 

ri= y dt, J dt2{u^{t,)u^{t2))e'^'2^^-)''- 
t t 
2-K (u^u_L)f^(t) At 
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t+At t+At 

t ti 



« -At / {u^ (r) (0)) e-*"(*)^dT (22) 



n, =- 



t+At ti 



t t 

oo 



-Ai y (u+ (r) u_L (0)) e-'"(*)^dT (23) 
^ I 



In this calculation we used the fastness of the noise 
{rn ^ At, r„ <C V^lj to substitute the integral in the 

exponent by O {t) {ti ~ and to extend the integration 
over the difference ti — 12 to infinite limits. The graph 6a 
connects Ni [t + At) to A'2 (t) , two others graphs connect 
A^i (t + At) to A^i (t) . Different signs in the contributions 
(PTjl and ([^ [221) are associated with the fact that the 
vertex at the upper line contains a factor — j, whereas it 
acquires the factor +i at the lower line. 



Collecting all contributions together, we arrive at the Master equations 



We rewrite themin a unified form: 



dNi 
~df 



n=n(t) 



(24) 



where 9 (x) is the step function equal to 1 at positive and at negative arguments. The master equation looks simpler 
when the variable Sz = — _ i ig used: 



dsz 
~dt 



2-K 



|f2 



-signal X 
2 ^ 



■ \n\ 



n=n{t) 



(25) 



The Fermi Golden Rule or the first Born approximation 
at a fixed moment of time can be applied since the fre- 
quency variation during the noise correlation is small and 
the perturbation caused by the noise in the correspond- 
ing stationary problem is weak. The perturbation the- 
ory is valid if <C |fi(t)|. Within the accumulation 
time interval the instantaneous frequency is of the same 
order of magnitude as the spectral width of the noise: 



\nit)\ 



The same inequality ensures that the fre- 



quency exceeds the width of the levels and the change of 
frequency due to the interaction with the noise. 

Now the question about the influence of the longitudi- 
nal noise onto the master equation is in order. First we 
demonstrate that correlations of the type (t)it|| (t')) 
do not change the master equation. Indeed, let us con- 
sider the influence of the longitudinal noise onto the dif- 
ference between Na (t -|- At) and Na (t) . In analogy with 
the case of the transverse noise, the contribution of one- 
longitudinal phonon line inside the interval (t, t -I- At) 
must be taken in account. This contribution does not 
depend on preceding evolution of the system. Therefore, 
it is the same as it would be in the absence of the trans- 
verse noise. Since the longitudinal noise does not produce 
transitions in the absence of the transverse noise, the to- 



tal contribution of 3 graphs of Fig. 6 for longitudinal 
noise is zero. 

V. RENORMALIZATION OF THE LZ GAP 



ti 



tl 12 



t2 tl 



t2 ti 



FIG. 7: Graphs containing mixed noise correlator (dash-dot 
line) and responsible for the LZ gap renormalization. 



The problem of mixed correlations between transverse 
and longitudinal noise is more subtle. The line of mixed 
correlation starts at one state, to say 1, and ends at an- 
other one (2) as shown in Fig. 7. The self-energy part 
associated with the graphs of the Fig. 7 reads: 
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(2 e 



(t') (0)) dt' ~ / {u^ (t') U|| (0)) dt' 



herm. conj . 



(26) 



This operator has the same form as the operator Aax 
in the interaction representation. Thus, at a time scale 
much longer than r„ , the mixed correlation renormalizes 
the LZ gap to a value 



A^A + i / {[u^ {0)])dt' 



(27) 



For the transformation of the integrals in equa- 
tion ([26|) into integral in equation (I27p we employed 
the time-translation invariance: (wy (— i) m_l (0)) = 
(0) wj^ (t)). Thus, the statement that one can ne- 
glect the action of the transverse noise within the LZ 
time interval is not completely correct: it is legitimate to 
neglect transverse-transverse correlations, but the mixed 
correlations can significantly change the LZ gap up to 
turning it into zero and changing its sign. 

The commutator entering equation ((27|) does not de- 
pend on the phonon occupation numbers, i.e. on tem- 
perature. It is instructive to express the renormalization 
of the LZ gap in terms of the phonon model (Section HIl 
equations (0 H El)): 



A 



V ^ 



.911 (q) .9-L (q) 



(28) 



VI. LONGITUDINAL NOISE 



Next we consider the action of the purely longitudi- 
nal noise within the LZ time interval. For this prob- 
lem the transverse noise will be ignored. If the noise 
is classical, the proper diagonal Hamiltonian is = 
^_ S7£t) ^ ,^||^ u ^ _|_ whereas the non-diagonal part is 

V — Ada;. In the interaction representation the non- 
diagonal part acquires the following form: 



y/(<) = a(|1) (2|e''''*o(^^- 



1u\\ \dr 



herm. conj. ) (29) 



Calculation of the transition probability is very similar to 
that considered above (see equation [20l and Fig. 2), but 
the vertexes correspond to Vi {t) given by ([29]) (we will 
call them A— vertexes) and instead of connecting pairs 
of noise amplitudes it is necessary to calculate average 

of a product exp ±2i u\^\^ {t j) dr j . Number of the 

j 

signs — in the exponent is equal to the number of signs +. 
Therefore, the dependence on the initial moment of time 
to vanishes. The Gaussian statistics allows to calculate 
the average of the product: 



J|exp 



±2i 



exp[— 2 



2fe-l 



*i 



dri. 



to 



to 



dT2kY^ (±(u||(ji)u|| (is)) (U||(j3)u|| (j4)) + •••)] 



(30) 



where summation is performed over all possible divisions 
of arguments tj into pairs. Each correlator vanishes if the 
modulus of corresponding time difference |r — r'| exceeds 
r„. The substantial range of integration over remaining 
variable {t + t') /2 is about tlz ~ A/il. Therefore the 
order of magnitude of the number obtained in the expo- 
nent ((5U1) after integration is ~ (^(u'^^'^ t„tlz^ ■ This 
value must be of the order or larger than 1 to ensure a 
significant change of the transition probability by the lon- 
gitudinal noise. This requirement is equivalent to the in- 

equahty ^"|^ ^ ^/ (At„) » tl. For the transverse noise 
the analogues criterion is much softer^^: (u^) > tl. 



The quantum longitudinal noise does not commute 
with itself at different moments of time. Therefore, it 
must be included into the non-diagonal Hamiltonian. 
Then, besides of A— vertexes one should consider the 
noise vertexes. We obtain the same estimate considering 
the contribution of an elementary graph with one noise 
line. However, for the fast longitudinal noise it is possible 
to find not only an estimate, but equations for the tran- 
sition amplitudes. First we prove that the contribution 
of the irreducible graphs^"*, which do not connect differ- 
ent branches of the Keldysh contour, is zero. Indeed the 
remote ends of such graphs are separated by the time in- 
terval of the order of t„ <C t^^- Therefore it is possible 
to integrate over the time difference at a fixed " center of 
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time" or "slow time" as we did in the case of the trans- 
verse noise. In contrast to the latter case, for longitudi- 
nal noise the integral does not depend on slow time since 
the longitudinal noise vertex in the interaction represen- 
tation docs not contain time-dependent phase factor (it 
connects identical states of the 2-state system). By the 
same reason it is the same for states 1 and 2. Therefore 
in a slow time scale much longer than r„ (but much less 
than tlz) the contribution of such graphs is proportional 
to unit operator for the 2-state system and can be com- 
pletely ignored. The irreducible noise graphs connecting 
different branches of Keldysh contour form a 4-pole ver- 
tex (see Fig. 8), which by the same reason does not de- 
pend on slow time and connects identical states at each 
branch of the Keldysh contour. We denote this vertex F. 



Po Po r P 

+ 



FIG. 8: Graphs containing the longitudinal noise only. Tri- 
angles correspond to the LZ gap A. 

a) , graphs with one noise loop. In the slow time scale they 
are equivalent to addition of a constant energy. 

b) , graphs with a line connecting Keldysh branches. 

c) , general graphic equation for P. See explanation in the 
text. 



Since the number of A— vertexes between vertexes F is 
arbitrary we need to extend the number of amplitudes in 
consideration. Namely, we define transition amplitudes 
Paf3,aoPg {t, to) as the average value of the operator |a) {f3\ 
at the moment t if the density matrix for 2-state system 
at the moment to was \ao) (/3ol- In the accepted approxi- 
mation these amplitudes obey a system of linear integral 
equations: 



t 



(31) 



Here P^fi aoPo (^'^o) denotes the transition amplitude in 
the absence of noise. The validity of equations (Plj) is 
limited by moderately strong noise. Otherwise the am- 
plitudes Pa0,ao0a (^j ^o) vary significantly at a time scale 
Tn- In this case more complicated equations with a non- 
local in time kernel F and amplitudes depending on 4 
time arguments must be used. 



VII. SOLUTION OF THE MASTER EQUATION 
AND NOISE DIAGNOSTIC 



The master equation (|25p allows an explicit solution: 



Sz {t) =Sz (to) exp - / / (r) dr 



(32) 



where / {t) = 27r (^(W)|f^(t)| + {n'^'n) 
-TTsignf] {t) ((W)|n(t)| - (^^^)-|a(t)| 



o(t) 



and g (t) = 
For classical 

noise g (t) — 0, and equation ((5^ reproduces the re- 
sult obtained in the reference (— ). At zero temperature 

W'^)-\n{t)\ = ^^^5 (*) = -^signf^ {t) f {t). In these 
two cases the measurement of occupation numbers or 
Sz (t) gives direct information on spectral power of noise. 
In classical case we find '^'^ {vV^)^fi(^t'f^ — ~^ln|sz|; in 
purely quantum case (T = 0) the relationship is slightly 



signO(i) 



more complicated: 27r (7777^)1^^^^! 

In general case it is possible to find both spectral func- 
tions / (t) and g (t) by performing two series of mea- 
surements with different initial states. Thus, the 2-state 
system is an ideal noise analyzer. 

Next we consider regimes of very fast and very slow 
(adiabatic) frequency sweeping. In the regime of fast 
sweeping (?/?/^) ^ tl the perturbation theory for equa- 
tion (|32p is valid. Indeed the integral f (r) dr can be 
rewritten in terms of spectral power as follows: 



o(t) 



to 



O(to) 



(33) 

The integral in the r.-h. side of equation ([55]) reaches its 
maximum value, equal to the average square fluctuation 
(7777^ -I- rjh'i'j at n (to) — —00 and {t) — +00. If the 
condition of fast sweeping is satisfied, the exponent in 
equation ([5^ can be expanded into a series over small 
noise parameter 7„ = (7777^ -f- r;^77) /(}. The variation 
Asz (t) = Sz {t) — Sz (to) is small at any time. In the 
leading approximation it reads: 

n(t) 

signw ^ 

o(to) (34) 



Aszit) ^ ^27rn-' I [{r^r^^)^^^isz{to)+^-Y^) 



In the opposite regime of slow (adiabatic) sweeping the 
noise parameter 7„ is large. In this case the exponents 
in equation (|32p vary very rapidly allowing asymptotic 
calculation of Sz- However, in adiabatic regime it is sim- 
pler to start directly with the Master equation . Ne- 
glecting the time derivative in it, we find the adiabatic 
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solution: 



g (t) _ signn (t) (wO|o(t)| - W'l)-m)\ 



fit) 



|0(t)| 



-m)\ 

(35) 

If the photon bath is in equihbrium with temperature T, 
equation (pSj) implies Sz (t) = — tanh . As it could 
be expected, at slow sweeping the two-state system adi- 
abatically accepts the equilibrium population with the 
temperature of the bath. This conclusion shows that in 
the case of the quantum noise one must be more care- 
ful with the asymptotic behavior of the time-dependent 
frequency than in genuine LZ problem or even in the 
analogues problem with the classical noise. In the lat- 
ter problems the linear approximation for CI (t) — fit 
was satisfactory. However, this approximation may be 
invalid for the quantum noise if the sweeping stops be- 
fore the frequency O (t) reaches the spectral width of the 
noise . In the opposite case the value Sz (t) saturates af- 
ter t — Tacc- In classical adiabatic case Sz (t) becomes 



zero after a short time Ttr = [\\Vj_\ / ^" j ■ similar 
time scale for the longitudinal noise was introduced by 
Kayanuma and Nakamura^*^. 

At the edge of the adiabatic regime 7„ ^ 1 the fast 
noise is still moderately strong, i.e. \/{r]T]^ <C r^^. It 
means that, when the noise becomes strong \/{r]ri^ > 
T~^, the system is already in deeply adiabatic regime. If 
the phonon bath is in equilibrium, the 2-state system also 
is in equilibrium with the noise. This equilibrium state 
is established in a time-independent, but strongly non- 
linear system. The interaction of the two-level system 
characterized by the time-independent frequency fl with 
the strong noise renormalizes the frequency and creates 
a finite width for each level. The situa tion b ecomes sim- 
pler in the limit of very strong noise \/(jpf) ^ ^ ^■ 
In this case the initial energy difference Q between levels 
can be neglected. The two states become equivalent and 
theiroccupation numbers are equal (1/2), i.e. Sz = 0. 
The same result can be obtained from the fact that, as 
we already argued, the very strong noise must be classi- 
cal. Equation (|35p can be considered as an interpolation 
between weak and very strong noise. Therefore, it gives 
a reasonable description of intermediate regime. 



time interval \t\ ^ tlz and to ignore the LZ gap A be- 
yond this interval. In this section we match the LZ so- 
lution modified by longitudinal noise inside the LZ inter- 
val with the solution of the problem with the transverse 
noise and A = (see Section IVII[) beyond this inter- 
val. For this purpose we choose a time scale ti such that 
tlz ^ ti ^ Tacc and first consider the solution of 
the noise problem with A = at t = —ti. For simpli- 
fication we accept to — -co. Since ti ^ Tacc it can 
be replaced by in the solution ([5^ with high precision 
~ ti/Tacc- Thus, at the left edge of the interval {—ti,ti) 
we find: 



Sz {~ti) ~ s( ' = exp ^- y / (t) dT^ Sz (-oo) 



(36) 



This value can be treated as an initial condition si ^ at 
t = — oo for the LZ problem with the longitudinal noise. 
If the solution of this problem is known, the value si"*^^ 
at t = +00 can be calculated. The information necessary 
to make this calculation effective is the knowledge of two 
numbers if there is no coherence in the initial system. 
The density matrix p(+^ at t = +oo is obviously a linear 
function of the initial density matrix p^~\ There exists 
a linear 4x4 matrix A performing this transformation: 



0^+' - A « o(-) 



(37) 



The requirement that Trp'^-' = 1 if Trp' ' = 1 implies 

the following equation: Aaa,fiu = 5^v If P12' — ^21' — 0, 
the equation l|37p results in following relationship be- 



tween si^^ and si ^ 



4+) = (Ai + A2) + (Ai - A2) 



(38) 



where we introduced abbreviations Ai and A2 for An^n 
and A22,22, respectively. If the longitudinal noise is ab- 
sent or sufficiently weak, the LZ values for A^ are: 



Ai = A2 = exp(-27r7i2) 



(39) 



VIII. TRANSITIONS IN THE PRESENCE OF 
THE LZ GAP AND NOISE 

As we demonstrated earlier, for the fast moderately 
strong noise, the effective time of the LZ transition due 
to the regular LZ gap tlz = A/f2 is much less than 

the accumulation time Tacc = {^Tri^ . Therefore it 
is possible to ignore the transverse noise within the LZ 



If {y^i ^ ^/ (At„), the longitudinal noise is weak 
enough to neglect the longitudinal-longitudinal correla- 
tions. Still the correlation of the longitudinal and trans- 
verse noise can significantly change the effective LZ gap 
(see equation ([27)) ) . 

The value si^^ from equation ([59]) serves in turn as 
initial condition at t = +0 for the Master equation. Its 
solution (I32I1 at t — +00 leads to the final result: 
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s,(+oo) = (Ai + A2)e-2-7^s,(-oo) + (Ai - A2)e-"T^^x 



Here 7^ = (u\) /tl. We remind that the occupation 
numbers are related to Sz as Ni^2 = 1/2 ± s^. Below 
we write the survival probability for the case when the 
longitudinal-longitudinal correlations can be neglected 



/f2 ^ (Ar„) In this case (—00) = 1/2, the 

values Ai^2 are determined by equation (|39|) and from 
equation ([3D]) we find: 



J 



Pi^i = 2 + ^'^''''^ (2e-2^'^^^ - 1) 



o 



(41) 



i]dn 



r 



To take in account the correlation between the longitu- 
dinal and transverse noise one should replace the LZ gap 
A in the expression for = A^/i7 by the rcnormalized 
value from equations p7|28p . 



IX. DISCUSSION AND CONCLUSIONS 

In the case of weak transverse noise or very fast sweep- 
ing 7j^ <C 1 equations ([Ml HT|) turn into the result ([55]) 
and the LZ survival probability, respectively. In the op- 
posite case of the strong transverse noise or slow sweep- 
ing 7j^ 3> 1 the occupation numbers accept their sta- 
tionary values at fixed instantaneous frequency indepen- 
dently on the value of LZ parameter 7^^. The classical 
noise corresponds to large phonon occupation numbers 



In this case the operators 77 and rf^ commute; all 
terms containing commutators (?7'7^)j^ — ('7^^)-fi '^^'^ ^e 



i,q. 



neglected. Then theory reproduces the result for classical 
fast noiseii. It is instructive to compare equation ([4T|) 
with the exact survival probability for T — Q obtained 
in the recent article by Wubs et al}^ . At zero tempera- 
ture the average value (r/'l' [t) rj {t')') as well as its Fourier 
transform turns into zero. This fact allows to calculate 
the integrals in equation ([^T|) . More physically visible 
way of obtaining the same result is to keep in mind that 
there is no live phonon at T = andonly spontaneous 
emission of phonons is possible. Therefore, if initially 
only the lower state was populated, the phonon cannot 
be emitted before the levels crossing. This considera- 



tion immediately gives si ■* = Sz (-co) = | 
i (2e~^^J^2 _ 1^. Employing general equation 
find the value (+00) — 5 [2 exp {—2tt (7^2 -|- 7_l)) 
and the survival probability: 



and 



(+) 



we 
-1] 



Pi^i = exp (-27r {j^z + 7±)) = exp[- 



27r 



1 



q 



511 (q) g± (q) . 



r + U {0)v{ (0)))] 



(42) 



r 



This result with precision of notations coincides with the 
exact result by Wubs et al,^^, equations (6-8), obtained 
without any limitations to the strength of noise and ra- 
tios of characteristic time scales. Surprisingly the multi- 
phonon processes as well as the longitudinal-longitudinal 



noise correlations do not contribute at all to the survival 
and transition probabilities even at very high noise in- 
tensity. At zero temperature such a high noise level can 
be reached only by enhancement of the coupling ampli- 
tudes. Though large coupling amplitudes are physically 
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implausible, as a mathematical model they are absolutely 
legitimate. The fact that these high-intensity processes 
do not play role in the transitions supports our specu- 
lations about possible extension of the master equation 
beyond their range of validity at least as a reasonable 
interpolation procedure. 

Our theory is relevant to molecular magnets, first of 
all because the condition of the noise fastness is per- 
fectly satisfied in the experiment. Indeed, the highest 
magnetic field rate used in the experiments with Fcg 
and Mni2 was 10^ Gs/s^r^. This rate corresponds to 
(l = 10^°s~^. The lowest temperature used in the cited 
measurements was about 0.05K. The dimensionless ra- 
tio of the value VTl to the smaller of the noise spectral 

widths is h\/ri/T ^ 10"°^. The estimates of the noise 
intensity is not so simple. The transverse noise is very 
weak for transitions with a large change of the spin pro- 
jection, for example, from +10 to -10, since the stan- 
dard coupling of the deformations to spin can change the 
spin projection only by ±1,±2. Therefore, the trans- 
verse noise for such a transition appears only in the 10- 
th order of the perturbation caused by acoustic phonons. 
However, for transitions with the small change of the 
spin projection the transverse noise is of the same or- 
der of magnitude as longitudinal one. The latter can be 
estimated from the known value of magnetic anisotropy 
energy Ea — QQK (for Feg) and the statistical weight of 
the phonons with the frequency smaller than anisotropy 
energy '--^ (Ea/En)^, where Ed is the Debye energy. An 
additional factor Ea/En comes from the square of cou- 
pling function |(7q|^, which is proportional to q at small 

wave vectors. Thus, ("f ) Ea{Ea/ Ed)"^ ^ Ea ~ 

hujg. This inequality shows that the condition of mod- 
erately strong noise is well satisfied, whereas the value 

/fi^^ is large. This fact together with the experi- 
mental fact that the hysteresis curve significantly narrows 
at T ^ 2 — 3K shows that the noise level is overestimated 
by this simple formula and more accurate theory is nec- 
essary. However, there is no doubt, that the thermal 



noise becomes important at a temperature of few Kelvin 
and that the noise is fast and quantum. Though the 
longitudinal noise does not produce transitions between 
diabatic levels, it is effective for transitions between adi- 
abatic levels^^. 

In conclusion, we derived the master equations for 
moderately strong fast quantum transverse noise. We 
obtained its solution at any moment of time. We demon- 
strated that the action of the regular LZ transition ma- 
trix element and the longitudinal-longitudinal noise cor- 
relation is limited by the LZ time scale, whereas the 
action of the transverse noise is accumulated during 
much longer accumulation time. We showed that the 
mixed longitudinal-transverse noise correlation leads to 
the renormalization of the LZ gap A. In the limiting 
case of adiabatic transverse noise the 2-state system adi- 
abatically follows its stationary state at an instantaneous 
value of frequency independently on the value of LZ pa- 
rameter. The separation of time scales allows to derive 
exact transition probability with the LZ gap, longitu- 
dinal and transverse noise taken in account simultane- 
ously. The transition probabilities depend explicitly on 
the noise commutator refiecting the quantum nature of 
the noise. It plays an important role especially in the 
adiabatic regime. In the extreme quantum regime at 
zero temperature our result coincides with exact result 
by Wubs et al .1^'. We argued that the strong- noise effects 
such as multiphonon processes and change of frequency 
appear only in the adiabatic regime for the fast noise and 
do not change substantially the transition probability. 
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